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■ Abstract. The pairwise entanglement, measured by concurrence and geometric phase in high dimensional 

free-fermion lattice systems have been studied in this paper. When the system stays at the ground state, 
their derivatives with the external parameter show the singularity closed to the phase transition points, and 
can be used to detect the phase transition in this model. Furthermore our studies show for the free-fermion 
model that both concurrence and geometric phase show the intimate connection with the correlation 
functions. The possible connection between concurrence and geometric phase has been also discussed. 

PACS. 03.65.Vf Phases: geometric; dynamic or topological; 03.65.Ud Entanglement and quantum nonlo- 
cality; - Q5.70.Fh Phase transitions: general studies 
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1 introduction 

The understanding of quantum many-body effects based 
on the fundamentals of quantum mechanics, has been rais- 
ing greatly because of the rapid development in quantum 
information theory [T]. Encouraged by the suggestion of 
Preskill p], the connection between the quantum entan- 
glement and quantum phase transition has been demon- 
strated first in ID spin-1/2 XY mode [3], and then was 
extended to more other spin-chain systems and fermion 
systems (see Ref [Ij for a review). Furthermore the de- 
coherence of a simple quantum systems coupled with the 
quantum critical environment has been shown the signifi- 
cant features closed to the critical points [HIS] . Regarding 
I these findings, the fidelity between the states across the 
transition point has also been introduced to mark the hap- 
pening of the phase transitions [B] . These intricate connec- 
tions between quantum entanglement and phase transition 
in many-body systems have sponsored great effort devoted 
to the understanding of many-body effects from quantum 
information point [T]. In general quantum entanglement 
as a special correlation, is believed to play an essential 
role for the many-body effects since it is well accepted 
that the non-trivial correlation is at the root of many- 
body effects. Although the ambiguity exists [7], quantum 
entanglement provides us a brand-new perspective into 
quantum many-body effects. However the exact physical 
meaning of quantum entanglement in many body systems 
remains unclear [8] . Although the entanglement witnesses 
has been constructed in some many-body systems [9], a 
general and physical understanding of quantum entangle- 
ment in many-body systems is still absent. 

On the other hand, the geometric phase, which was 
first studied systemically by Berry [10] and had been re- 



searched extensively in the past 20 years [TT], recently 
has also been shown the intimate connection to quantum 
phase transitions [T21[T31[Ill[Ill[ini[lL,18 (or see a recent 
review Ref. [19]). This general relation roots at the topo- 
logical property of the geometric phase, which depicts the 
curvature of the Hilbert space, and especially has direct 
relation to the property of the degeneracy in quantum sys- 
tems. The degeneracy in the many-body systems is critical 
in our understanding of the quantum phase transition . 
Thus the geometric phase is another powerful tool for de- 
tecting the quantum phase transitions. Moreover recently 
geometric phase has been utilized to distinguish differ- 
ent topological phases in quantum Hall systems [21\, in 
which the traditional phase transition theory based on the 
symmetry-broken theory is not in function [22| . 

Hence it is very interesting to discuss the possible con- 
nection between entanglement and geometric phase, since 
both issues show the similar sensitivity to the happening 
of quantum phase transition. Recently the connection be- 
tween the entanglement entropy and geometric phase has 
first been discussed with a special model in strongly cor- 
related systems; the geometric phase induced by the twist 
operator imposed on the filled Fermi sphere, was shown to 
present a lower bound for the entanglement entropy |23j . 
This interesting result implies the important relation be- 
tween quantum entanglement and geometric phase, and 
provides an possible understanding of entanglement from 
the topological structure of the systems. In another way 
the two-particle entanglement was also important f3]. Es- 
pecially in spin-chain systems two-particle entanglement 
is more popular and general because of the interaction 
between spins, and furthermore the quantum information 
transferring based on spin systems are generally depen- 
dent on the entanglement between two particles [24^. So 
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it is a tempting issue to extend this discussion to the uni- 
versal two-particle entanglement situation. 

For this purpose the pairwise entanglement and geo- 
metric phase are studied systemically in this paper. Our 
discussion focuses on nearest-neighbor entanglement in 
the ground state in free-Fermion lattice systems because 
of the availability of the exact results. By our own knowl- 
edge, this paper first presents the exact results of entangle- 
ment and geometric phase in higher dimensional systems. 
In Secl2]the model will be provided, and the entanglement 
measured by Wootter's concurrence is calculated by intro- 
ducing pseudospin operators. Furthermore the geometric 
phase is obtained by imposing a globe rotation, and its 
relation with concurrence are also discussed generally. In 
SecEl we discussed respectively the concurrence and geo- 
metric phase in 2D and 3D cases. Finally, the conclusion 
is presented in Sec|4l 

2 Model 

The Hamiltonian for spinlcss fermions in lattice systems 
reads 

^ = E ^.^^y^j + 2^ABuc] + h.c), (1) 

ij 

in which is fermion annihilation(creation) operator 
and L is the total number of lattice sites. The hermi- 
tity of H imposes that matrix A is Hermit and B is an 
anti-symmetry matrix. The configuration of lattice does 
not matter for Eq. ^ since our discussion focuses on the 
general case and available exact results. This model obvi- 
ously is solvable exactly and can be transformed into the 
free Bogoliubov fermionic model. So it is also called free- 
fermion model. By Jordan- Wigner transformation |27| one 
can convert the spin-chain systems into spinless fermions 
systems, in which the physical properties can be readily 
determined. Therefore an alternative approach is neces- 
sary by which one can treat solvable fermion systems of 
arbitrary size. The model Eq. ([T]) serves this purpose. 

Without the loss of generality we assume A and B to 
be real [55] . An important property of Eq. ([1]) is 

L 

[i/,n(l-2cK)] = 0. (2) 

i 

This symmetry would greatly simplify the consequent cal- 
culation of the reduced density matrix for two fermions. 
One can diagonalize Eq. ([T]) by introducing linear trans- 
formation with real gti and h\a |25| 

1 ^ 

■>lk = ^j=^9\LiCi + h]^ic\, (3) 

in which the normalization factor \/y/L have been in- 
cluded to ensure the convergency under the thermody- 
namic limit. After some algebra, the Hamiltonian Eq. ([1]) 
becomes 

= E ^k??l77k + const. (4) 

k 



in which is the common eigenvalue of the matrices 
{A - B){A + B) and [A B){A - B) with the corre- 
sponding eigenvectors (/)ki = 5ki + ^ki and V'ki = 5ki — ^ki 
respectively (see Ref. [IS] for details). The ground state is 
defined as \g), which satisfies the relation 

^/klff) = (5) 
With respect to fermi operator 77k, one has relations 

E 3ki5k'i + /iki^k'i = <5^^k 
i 

E ffki/ik'i + ^ki5k'i = (6) 
i 

Furthermore the requirement that {(/)fe,V/c} and {V'fc,VA;} 
be normalized and complete, reinforce the relations [25j 

E 5ki5kj + ^ki/lkj = 5ii 
k 

E 5ki^kj + /lki.9kj = (7) 
k 

With the help of these formula above, one obtains 

Ci = E 5ki'/k + /iki?7^ (8) 

which would benefit our calculation for the correlation 
functions. 



2.1 Concurrence 

The concurrence, first introduced by Wootters [26] for the 
measure of two-qubit entanglement, is defined as 

c = max{0, Ai — A2 — A3 — A4}, (9) 

in which = 1, 2, 3, 4) are the square roots of eigenval- 
ues of matrix R — p{ay (g) a'^)p{ay (g) tr^) with decreasing 
order. Then the critical step is to determine the two-body 
reduced density operator p. The reduced density operator 
Pij for two spin-half particles labeled i,j can be written 
generally as, 

I 4 

Pij = ti-ijP = 4 E Po'-.P^i ® '^f ' (10) 

a, 13=0 

in which p is the density matrix for the whole system and 
cr° is the 2x2 unity matrix and i7"(a = 1,2,3) are the 
Pauli operators ct^,(T^,(t^, which also the generators of 
SU{2) group. Pq,/3 = tr[(T"tT?/9ij] = {^f^j) is the correla- 
tion function. With the symmetry Eq. Q, one can verify 
that only Poo, P03,P30,Pii,P22,P33,Pi2,P2i are not vanish- 
ing. After some efforts, one obtain 

c = max{0, c/, c//}, (11) 
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in which 



c/ 



V (1 +P33)^ - (P30 +P03)^] 



C/7 



-P22I - V(l - -P33)^ - (P30 - P03)^]-(12) 



In order to obtain the reduced density operator for 
two fcrmions, it is crucial to construct SU{2) algebra for 
the ferniions in lattice systems. In ID case, the Jordan- 
Wigner (JW) transformation is available [23 US [13, 25 . 
For higher dimension cases the JW-like transformation 
has been constructed by different methods [28]. However 
the transformation is very complex and the calculation is 
difficult. Hence instead of a general calculation, we focus 
on the nearest neighbor two lattices in this paper. In this 
situation, the SU (2) algebra can be readily constructed 



ometric phase(GP). Then we have Hamiltonian with pa- 
rameter (f) 

H{<j>) = ^4^,,cj + -(^Bycje^.^ + (15) 



and the ground state becomes \g{(p)) — R{(j))\g). GP is 
defined as [TO] 



73 = -*/ d^igm-^im 



(16) 



Regarding to Eg . (fT5| . one only require t/i = tt for a cycle 
evolution. Hence one has 7ff = f Si J2k ^ki = r Si 7si- 



'i+l 



^i+l 



^i+i 



.- = (ar-zaf)/2 = ci 
a? = 2c\ci - 1 
_i)/2 = (2ctci - l)4i 
_,)/2 = {2clc, - l)ci+i 



^i+i 



= 2C: 



t 

i+1 



Ci+1 



- 1 



(13) 



in which i + 1 denotes the nearest neighbor lattice for site 
i. This point can be explained as the following. The diffi- 
culty for the JW transformation in higher dimension case 
comes from the absence of a natural ordering of particles. 
However when one focuses on the nearest neighbored par- 
ticle, this difficulty does not appear since for a definite 
direction the nearest neighbor particle is unique (for non- 
nearest neighbored case one have to consider the effect 
from the other particles). Then the correlation functions 
for the ground state are in this case 



Poo 



1,P30 = 1 - 



Pii = 



P22 



P33 



k 



L ^ 

k 



hli, Pas = 1 - 



L ^ 

k 



'-k(i+l) 



k(i+l) 



.9k(i+i)); 



5ki)(^: 



k(i+i) 



.9k(i+l)j 



(i-T>:/^ki)(i 



L ^ 

k 



--j-2 'iki/jk(i+l)5k'i.gk'(i+l) - 'lki5ki/jk'(i+l)5k'(M 
k.k' 



2.3 GP vs. Concurrence 



At a glance of Eq. (fT^ and Eq. ([TB]) , GP and concurrence 
both are related directly to correlation functions. Hence it 
is tempting to find the relation between the two quantities, 
which would benefit to the understanding of the physical 
meaning of concurrence. 

According to Eqs. fT^ and ([H]), the following inequal- 
ity can be obtained (see Appendix for details of calcula- 
tions) 

1 



CI < J^ilgi + lg(.i+i)) ■ 



cii < 1 



7s(i+i)) 



V(l + P33)^ - {P30 + PQs)"^ 



2L^7 



For the first inequality, a much tighter bound is difficult 
to find. While if the average of c// over all site i is con- 
sidered, c// < 1 - Y^iilgi - 7g(i-n))^- Fortunately in 
the following examples c/ is always negative. Although the 
existence of this defect, in our own points, the relation be- 
tween GP and concurrence have been displayed genuinely 
from the inequality above. 



3 GP and Concurrence in Higher Dimensional 
XY model 

The previous section presents the general discussion of GP 
and concurrence in free fermion lattice system Eq.([T]). In 
this section a concrete model would be checked explicitly, 
of which the Hamiltonian is 



1) 



P12 =P21 =0 



(14) 



(18) 



2.2 Geometric Phase 

Following the method in Refs. [T^ll3j . one can introduce 
a globe rotation R{(j>) = exp[i0^,- C;Ci] to obtain the ge- 



in which denotes the nearest- neighbor lattice sites 
and Ci is fermion operator. This Hamiltonian, first intro- 
duced in Ref. [29j, depicts the hopping and pairing be- 
tween nearest-neighbor sites in hypercubic lattice systems, 
in which A is the chemical potential and 7 is the pairing 
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potential. Eq. (fT8|) could be considered as a d-dimensional 
generalization of ID XY model. However for d > 1 case, 
this model shows different phase features [29] . 

The Hamiltonian can be diagonalized by introducing 
the d-dimensional Fourier transformation with periodic 
boundary condition in momentum space [29j 



under the thermodynamic limit 



iZ\k(4c-k - h.c), 



(19) 



in which = Sa=i '-'^^ — X and Z\k — 7 X]q=i si'^ ^a- 
With the help of Bogoliubov transformation, one obtains 



= ^ 2ylk77^^k + const. 



(20) 



in which ylk = V^k^ ^k- Based on the degeneracy of 
the eigenenergy ylk = 0, the phase diagram can be de- 
termined clearly [29]; When d — 2, the phases diagram 
should be identified as two different situations; for 7 = 0, 
the degeneracy of the ground state occurs when A G [0, 2], 
whereas the gap above the ground state is non-vanishing 
for A > 2. However for 7 7^ three different phases can 
be identified as A = 0, A e (0, 2] and A > 2. The first two 
phases correspond to case that the energy gap above the 
ground state vanishes, whereas not for A > 2. One should 
note that A = means a well-defined Fermi surface with 
kx — ky ± IT, whose symmetry is lowered by the pres- 
ence of A term. For d — 3 two phases can be identified as 
A € [0, 3] with the vanishing energy gap above the ground 
state and A > 3 with a non-vanishing energy gap above 
ground state. In a word the critical points can be iden- 
tified as Ac — d{d = 1,2,3) for any anisotropy of 7, and 
A = for c? = 2 with 7^0. One should note that since 
the 7^ dependence of ylk, the sign of 7 does not matter. 
Hence the plots below are only for positive 7. 

The correlation functions between nearest-neighbor lat- 
tice sites would play a dominant role in the transition 
between different phases because of the nearest-neighbor 
interaction, similar to the case in XY model [3]. Then it is 
expected that the pairwise entanglement is significant in 
this model. In the following, concurrence for the nearest- 
neighbor sites of ground state is calculated for d = 2,3 
respectively. The geometric phase of ground state is also 
calculated by imposing a globe rotation R{(t>). our calcu- 
lation shows that both quantities show interesting singu- 
larity closed to the boundary of different phases. 



3.1 Concurrence 



For d > 1 case, the nearest-neighbor lattice sites appear 
in different directions. In order to eliminate the depen- 
dence of orientations, the calculation of correlation func- 
tions Eqs. (|14p is implemented by averaging in all direc- 
tions. With the transformation Eq. (|13p . one can determine 



Pll 



y / ]Jdfca(Z\fe ^ Sinfca - tfc ^ COsA;a)/ylfc 
^ «7r d d d 

P22 = - w / Y\dkaiAk^ sinka +tk^ COS ka)/Ak 



a=l 



P12 ^ P21 = 

P03 = P30 = Pa = 



(2 



2 ,pii +P22.2 , fPii -P22.2 

P33=P3-( ^ ) +( 7, ) 



(21) 



d = 2 Our calculation shows that c/ is negative. So in 
Fig. (TJ only c// and its derivative with A are numerically 
illustrated. In order to avoid the ambiguity because of the 
cutoff in the definition of concurrence, the derivative of 
cii with A is depicted in all region whether c// positive or 
not [7j . Obviously the singularity for dcu / dX can be found 
at the point A = 0, 2 respectively, which are consistent 
with our knowledge about phase transitions. 

d = 3 Similar to the case of d = 2, our calculation 
shows c/ < 0. Only c// and its derivative with A are nu- 
merically displayed in FiglS] Different from the case of 
d = 2, no singularity of the first derivative of c// with A 
is found at A = 3. While a cusp appears at A = 1. A fur- 
ther calculation demonstrates that the second derivative 
of c// is divergent genuinely at exact A = 3, as shown 
in FigsHI^c). which means the phase transition at this 
points. Furthermore our numerical calculations show that 
d^cii/dX^ is finite at A = 1, as shown in Figsl2l^b). Hence 
one cannot attribute this feature to the phase transition. 
The similar feature has been found in the previous stud- 
ies [SllTllS^. However the underlying physical reason is 
unclear in general. But this special feature is not unique 
for concurrence; van Hove singularity in solid state physics 
displays the similar feature, which is because of the van- 
ishing of the moment-gradient of the energy. Although 
we cannot established the direct relation between these 
two issues because of the bad definition of the moment- 
gradient of the energy when degeneracy happening, we 
affirm that this feature is not an accident and the under- 
lying physical reason is still to be found. 

In a word the discussion above first demonstrates the 
exact connection between concurrence and quantum phase 
transitions in high-dimcnsional many body systems. How- 
ever a question is still open; what the physical interpreta- 
tion of concurrence is in many-body systems. In this study, 
we includes the negative part of c// to identify the phase 
diagram in free-fermion systems. In general, it is believed 
that the negative c// means no entanglement between two 
particles and then include no any useful information about 
state. But from the discussion one can note that the omis- 
sion of the negative part of c// would lead to incorrect 
results. Moreover, for 7 = 0, our calculations show that 
c/, c// always are zero, and so one cannot obtain any the 
phase transition information from pair wise entanglement 



H.T. Cui, et.al.: Pairwise Entanglement and Geometric Phase in Higher Dimensional Free-Fermion Lattice systems 5 



in this case. Further discussions will be presented in the 
final part of this paper. 

3.2 Geometric Phase 

Geometric phase manifests the structure of Hilbert in the 
system and has intimate relation to the degeneracy. GP, 
defined in Eq. (fT6|) by imposing a globe rotation on 
ground state \g) is calculated in this section. After some 
algebra, one obtains 

d = 2 In FigO 7g and its derivative with A are dis- 
played explicitly. Obviously one notes that djg/dX shows 
the singularity closed to A = 0, 2, which are exactly the 
phase transition points of Hamiltonian Eq. (fT8|) . An inter- 
esting observation is that closed to these points, both GP 
and concurrence c// show the similar behaviors. 

d — S GP and its derivative are plotted explicitly in 
Fig. ([4]). One should note that there is a platform below 
A = 1 for djg/dX, as shown in Fig|4l^a), but a further cal- 
culation shows that d^jg/dX^ is continued (FiglH^b)) and 
d"fg/dX has no divergency at this point. This phenomena 
is very similar to the case of concurrence (see Figl2jb, c)). 
As expected, d'^^fg/dX^ is divergent at exact A = 3, which 
means a phase transition happens at this point. Together 
with respect of the case of d = 2, it makes us a suspect that 
GP and concurrence in our model have the same physical 
origination. 

Furthermore for 7 = 0, GP fails to mark the phase 
transition too. This is similar to the case of concurrence, 
but has different physical reason. The further discussion 
is presented in the next section. 



4 Discussion and Conclusions 

The pairwise entanglement and geometric phase for ground 
state in d-dimensional (d = 2,3) free-fermion lattice sys- 
tems are discussed in this paper. By imposing the trans- 
formation Ea. (jl3p . the reduce two-body density matrix 
for the nearest neighbor particles can be determined ex- 
actly for any dimension, and the concurrence is also calcu- 
lated explicitly. Furthermore geometric phase for ground 
state, obtained by introducing a globe rotation i?(0), has 
also been calculated. Given the known results for XY 
model [3l[T2l[13], our calculations show again that both 
GP and concurrence display intimate connection with the 
phase transitions. Moreover an inequality relation between 
concurrence and geometric phase is also presented in Eq. 
P7|) . The similar scaling behaviors at the transition point 
A = 3 has also been shown in Figs.[5l These facts strongly 
mean the intimate connection between the two items. This 
point can be understand by noting that both of them are 
connected to the correlation functions, as shown in Eqs. 
(HH) and (HH). 



An interesting point in our study is that in order to 
obtain all information of phase diagram in model Eq. (jl8p , 
the negative part of c// has to be included to avoiding 
the confusion because of the mathematical cutoff in the 
definition of concurrence [7] . In general, it is well accepted 
that the negative part of c// gives no any information of 
quantum pairwise entanglement, and then is considered to 
be meaningless. However, in our calculation, the negative 
part of c// appears as an indispensable consideration to 
obtain the correct phase information. This point means 
that the pairwise entanglement does not provide the all 
information about the system since the two-body reduced 
density operator throw away much information. 

As for the geometric phase, defined in Eq. (fTB)) . it is 
obvious that 7^ can tell us the happening of phase transi- 
tion at the point, where 7^ display some kinds of singular- 
ity. However it cannot distinguished the degenerate region 
from the nondegenerate, as shown in Figs. [3] and [H Re- 
cently GP imposing by the twist operator in many-body 
systems is introduced as an order parameter to distinguish 
the phases |T7lfT8| . For the free-fermion lattice system, this 
GP have also calculated and shows the intimate connec- 
tion with the vanishing of energy gap above the ground 
state. However the boundary between the two different 
phases becomes obscure with the increment of dimension- 
ality in that discussion [T7], and moreover it cannot dis- 
tinguish the phase transition not come from the degener- 
acy of the ground-state energy. While the geometric phase 
imposing by the globe rotation R{(j)) clearly demonstrate 
the existence of this kind of phase transition, as shown in 
Fig|31 whether originated from the degeneracy or not. In 
fact this point can be understood by noting the intimate 
relation between 7g and correlation functions. It maybe 
hint that one has to find different methods for different 
many-body systems to identify the phase diagram. 

Although the intimate relationship of concurrence and 
GP with phase transitions in the model Eg. p^ . a excep- 
tional happens when 7 = 0, in which c/, c// are zero and 
GP is a constant independent of A. From Eq. p8|) . 7 = 
means the hopping of particles is dominant, and the posi- 
tion of particle becomes meaningless. Since the calculation 
of concurrence depend on the relative position of lattice 
site, the pairwise entanglement is disappearing. However 
one could introduce the spatial entanglement to detect the 
phase transition in this case [31] . For GP, 7 = means the 
emergency of new symmetry. One can find c\ci,H] = 
in this case, which leads to the failure of R{4i) for construc- 
tion of nontrivial GP. 

Finally we try to transfer two viewpoints in this pa- 
per. One is that concurrence and geometric phase can be 
used to mark the phase transition in many-body systems 
since both of them are intimately connected to the cor- 
relation functions. The other is that concurrence and the 
geometric phase are connected directly by the inequal- 
ity Eq. P?)) . Then it is interesting to extend this relation 
to multipartite entanglement in the future works, which 
would be helpful to establish the physical understanding 
of entanglement. 
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APPENDIX 

For the first inequality, one sliould note 

bii + P22I 

k k 

< I Y.^hl, + = + 7^(i+i)). (23) 

k 

From inequality -^/a;^ — > |a;| — \y\{\x\ > \y\), one re- 
duces 

V(1+P33)^-(P30+P03)^] > |l+P33|-b30+P03|- (24) 

Then one obtains 

ci < ;^(7gi + 7s(i+i)) + ^(|P30 +P03I - ll+PssD- (25) 

However a much tighter bound is difficult to decide be- 
cause of the complexity of P33. 

For the second inequality, it can be obtained easily by 

observing 

P33<l-;^(7,V7|(i+i)), (26) 

in which we have used the relation 2ab < + b'^. Then 
1 — P33 is non- negative and 

^ S^\U I , P33 - 1 

c// = ^ 2^ \hki9k{i+i)\ + — ^ — 

< J $^(^ki + 5k(i+i)) - ^^(Tfli - 

k 

^ 1 + ;^(7si - 7fl(i+i)) - ^^(7si - la{i+i)p7) 
in which l/LJ2i^ 9ui = ^~ ^/^J2k ^ki used. 
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Fig. 1. cii (O) and its derivative with A (A) vs. A when 
d = 2. We have chosen 7 = 1 for this plot. 
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Fig. 3. 7g (O) and its derivative with A (A) vs. A when ( 
We have chosen 7 = 1 for this plot. 
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Fig. 2. cii (O) and its derivative with A (A) (a) vs. A when 
d = 3. We have chosen 7 = 1 for this plot. The second deriva- 
tive of Cii with A are also displayed in this plot and focus the 
points closed to A = 1 (b) and A = 3 (c). 
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Fig. 4. 7g (O) and its derivative with A (A) (a) vs. A when d = 
3. We have chosen 7 = 1 for this plot. The second derivative' 
of 7g with A are also displayed in this plot and focus the points 
closed to A = 1 (b) and A = 3 (c). 
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Fig. 5. The scaling of GP and concurrence for 3D case closed 
to the critical point Ac = 3. We have chosen 7 = 1 for this plot. 



